Sub-Rayleigh Quantum Imaging Using Single Photon Sources 
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We propose a technique capable of imaging a distinct physical object with sub-Rayleigh resolution 
in an ordinary far-held imaging setup using single-photon sources and linear optical tools only. We 
exemplify our method for the case of a rectangular aperture and two or four single-photon emitters 
obtaining a resolution enhanced by a factor of two or four, respectively. 
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In modern quantum optics, there is a great variety of 
proposals trying to improve different aspects of the image 
formation process, commonly summarized by the field of 
quantum imaging. Today, this fast growing field ranges 
from early Ghost imaging [l[ , sub-wavelengthphase mea- 
surements @, [H to quantum lithography [J, H, 0] , quan- 
tum microscopy 0, H S E3] and many more (see e.g. llj]). 
Though all proposals commonly aim to overcome the 
classical boundaries of image formation, only few im- 
prove the spatial resolution itself, i.e. the ability to image 
a physical object while overcoming the Rayleigh [l2| or 
Abbe limit [l3| of classical optics. 

So far, in quantum imaging sub-classical resolution 
has been achieved by using sources of entangled pho- 
tons [|| , but it was also shown recently that init ially 
uncorrelated light can be used for that purpose 1(J H|. 
All those methods exploit second (or Nth) order corre- 
lations between two (or N) photons, i.e. quantum in- 
terferences between two- (or N-) photon amplitudes, to 
surmount the classical boundaries. Hereby, it is yet a 
challenge to implement sub-Rayleigh quantum imaging 
using linear optical tools only. 

In this letter, we propose a method of imaging a physi- 
cal object, e.g. an aperture, beyond the classical Rayleigh 
resolution using linear optics. Our scheme involves N 
uncorrelated single photon emitters serving as a non- 
classical light source and N detectors performing correla- 
tion measurements. By placing the N detectors at differ- 
ent positions in the Fourier plane of the object we avoid 
the use of multiphoton absorption techniques. Moreover, 
using a lens in the Fourier plane our setup is also capable 
to reproduce the object in the image plane of the lens. 
We exemplify our method for the case of two single pho- 
ton emitters. By exploiting two-photon interferences we 
show that this scheme allows to image the object with 
sub-Rayleigh resolution, i.e., with a resolution enhanced 
by a factor of two with respect to the classical case. In the 
same way, sub-Rayleigh resolution enhanced by a factor 
of four is obtained for AT = 4 emitters and using N = 4 
detectors. By extending this scheme, we show that the 
same results are also obtained for different objects, e.g., 



in case of a grating with A" slits. 

Only recently, first experiments were able to observe 
higher order interferences of photons emitted by sin- 



gle trapped atoms [15|, [lg, [171, [18|. These observa- 



tions stand in a long line of experiments using single 
photon sources to investigate interference phenomena of 
single- and multi-photon amplitudes. After the early 
demonstration of first-oder interferences of light scat- 
tered by two ions [22| . two-photon interferences have 
been observed and quantified by now in several sys- 
tems 0, [H [H H, EI 0, H3- Generally, quantum 
interferences may appear if a measurement does not allow 
to reveal information about the particular path a quan- 
tum state has propagated (so called Welcher Weg infor- 
mation). For example, if we detect N (N > 1) photons 
by N different detectors placed in the far-field region of a 
source consisting of A~ single photon emitters, there are 
generally A! possibilities how the A" photons may prop- 
agate from the A atoms to the A" detectors [r| • Hereby, 
the N\ quantum paths may differ from each other by an 
optical phase, leading to destructive or constructive in- 
terferences between the A'-photon amplitudes. This can 
be fruitfully exploited, e.g., to obtain information about 
the spatial distribution of the source even if of dimen- 
sions smaller than the optical wavelength A 1,0. In 
the following, we will extend this concept by introducing 
a distinct physical object between source and detectors 
which we aim to image with sub-Rayleigh resolution in 
an ordinary far-field imaging setup. 

The configuration used in our imaging scheme is shown 
in Fig. [T] It consists of two identical two-level atoms 
at Ri and R2 initially excited by a single laser 7r pulse 
(where we denote the initially excited state by |ei,e2)). 
After the spontaneous decay of both atoms, the two pho- 
tons may be recorded by two detectors placed at ri and 
Y2 in the far-field region of the source. We next place 
a physical object between the atoms and the detection 
plane; in a first step we consider a rectangular aperture 
with opening height a and width b. In a single successful 
measurement cycle, the two photons emitted by the two 
atoms pass both by the aperture (distinct to [H) and are 
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FIG. 1: (color online) Setup used for quantum imaging: in a 
single measurement cycle both emitters at Ri and R2 each 
emit a single photon. Between the source and the detection 
plane, the field is diffracted by an object, e.g. a rectangular 
aperture of height a and width b. The two detectors are placed 
in the far-field region of the aperture at positions ri and r2. 



registered by the two detectors. For the sake of simplicity 
the two atoms and the two detectors are located in the 
x-y-plane, parallel to the plane of the aperture, where in 
particular R\ z = R% z := R z and r\ z = r 2i := r z . More- 
over, we consider coincident detection only 23]. Cor- 
relating the two detector events we measure the second 
order correlation function 

G( 2 )( ri ,r 2 ) = 
(e 1 ,e 2 |^(-)(r 1 ) J E(-)(r 2 )^(+)(r 2 )^(+)(r 1 )|e 1) e 2 ), (1) 

which, in the present setup, can be understood as the 
(unnormalized) joint probability to find one photon at ri 
and another one at r 2 . Hereby, £ (+) (r*) = [ £;( ~ ) ( r i)] t 
denotes the positive frequency part of the electric field 
amplitude at point r^. 

From classical diffraction optics (see e.g. [2^) we know 
that the electric field E{r) is diffracted at an aperture A 
placed between the source and the detection plane. The 
disturbed field U(r) can be calculated employing Fresnel- 
Kirchhoff diffraction theory. Using standard Fresnel 
approximations between source, aperture and detection 
plane we obtain for the diffracted field at the following 
expression (see (23 |) 
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where A denotes the initial amplitude of the electric 
field, k = y, ro is a vector in the plane of the aper- 
ture (see Fig. [T]) and is a vector consisting of the x- 
and y-components of Rj, ro and Tj, respectively (with 
7r = j, 0,i; j, i — 1,2). As one can see from Eq. ((2]), 
the problem is separated into the propagation from the 
source to the aperture (i.e. from Rj to ro) and further 
from the aperture to the detection plane (i.e. from ro to 



r;). Since we assume that in our setup the far- field con- 
dition is fulfilled we can limit ourselves to Fraunhofer- 
diffraction, i.e., 



n\ » y |po| 2 , 



(3) 



(though the Fresnel-integral @ can be solved numeri- 
cally without this approximation), fn this limit, one can 
carry out the integral over A in Eq. ([2]) to obtain the 
final expression for the disturbed field 



U(ri,R j ) = — e 
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In Fig. [TJ two atoms contribute to the electric field at 
Ti, each giving rise to a field U(ri, Rj) of the form given 
in Eq. ([4]). We can thus write the total positive frequency 
part of the field contributing to the correlation signal at 
r, as 



E(+\r i ) = ^=(U(v h R 1 )\g) 1 



e |+C/(r 4 ,R 2 )| ff ) 2 (e|), (5) 



where the atomic operator |g)j(e| describes the de- 
excitation of the jth atom. With i?( + )(rj) at hand, we 
can write the second order correlation function, Eq. |T]), 
in the form 



G^(r l7 r 2 ) = 
JWn.RiMra.Ra) 



(6) 



In classical optics, using a coherent light source and 
in the limit of Fraunhofer-diffraction, we know that a 
rectangular aperture with opening height a and width b 
gives rise to the following classical intensity diffraction 
pattern at point r in the far field (see [2J]) 
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where R z is the distance between the source and the aper- 
ture. From this expression we can easily recover Abbe's 
criterion for resolving the aperture from the intensity dis- 
tribution I(r) in the far-field, i.e. in the Fourier plane: 
this is the case only if it is possible to measure the in- 
tensity diffraction pattern in the range — 27r < 4^ < 2tt 
(a <-> b), i.e., if we obtain the first diffraction order. In 
the words of Abbe: we can reconstruct an image from 
an object if and only if the first diffraction order in the 
Fourier plane is at least visible [I3| ■ Otherwise, if we try 
to image an object of smaller sizes so that the diffrac- 
tion pattern in the Fourier plane increases beyond that 
range, the image will start to blur. Note that Abbe's 
criterion can be used to define the resolution limit in any 
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imaging technique (classical or non-classical) . The crite- 
rion enables us in the following to compare the classical 
limit of I(t) with the two-photon correlation signal of our 
proposed imaging scheme. 

Using Eqs. and ©, it is possible to explicitly cal- 
culate the second order correlation function G^ 2 ^(ri, r^). 
For the sake of simplicity, we will fix the detectors and the 
source in the x-z-plane, i.e., r iy = Rj y = 0. In difference 
to the classical intensity pattern, we have for (ri , r 2 ) 
to determine two parameters when performing a mea- 
surement, namely and r 2 . In addition, we have to fix 
the position of the two single photon emitters with re- 
spect to the aperture. The latter determines the phase 
shift between the different photon paths leading from ei- 
ther of the two emitters to the object. For example, if we 
choose Ro — R\ + n^r- and Ira I = ri := r, we obtain 

\ ( 8A 2 r 2 \ 2 . 9 kar 

G( )(ri ^H ^w ± (o ) - sm 2 ^' (8) 

where B+{r) = r 2 + ^f [5_(r) = |(r 2 -£^§)] holds for 
T2 X = + r i x [ r 2 x = -''lj,]- Comparing the modulation of 
the classical intensity pattern I(r) from Eq. (J7|) with the 
one obtained for the G^-function in Eq. ((8]), one can see 
that the latter oscillates twice as fast as in the classical 
case by keeping a contrast of 100%. Here, the additional 
degrees of freedom present in the G^-measurement en- 
able to select only those two-photon amplitudes which 
accumulate the phase |^ twice as fast as in the classical 
case. The increase of the modulation frequency by a fac- 
tor of two implies, according to the Abbe criterion, that 
sufficient information is available in the Fourier plane to 
reconstruct the aperture even if measuring only half of 
the range needed for the classical imaging technique. 

The results found in the limit of Fresnel and Fraun- 
hofer approximations for the case of N = 2 emitters using 
N = 2 detectors can be extended for the case of N > 2. 
For example, if we consider N = 4 single photon emitters 
located at positions Rx x — —w^-, i? 2;c = 0, i?3 x = \ % L , 
i?4 x = ttx*- and choose for the four detectors positions 
I raj = r, r 3m = -r + n^, r 4x = r + §fj, where again 
T\ x := r, we obtain for the fourth order correlation sig- 
nal [25| the following expression 

G( 4 W 2 ,r 3 ,r 4 )ocsin 2 4^. (9) 

2r z 

As before, we find that it is possible to image the object 
with sub-Rayleigh resolution. However, using four emit- 
ters and four detectors, the resolution is now enhanced 
by a factor of four with respect to the classical case. 

Note that our scheme also allows to reproduce the ob- 
ject in the image plane of a lens placed in the Fourier 
plane of the object. Fig. recalls the initial setup under 
investigation in case of N — 2. We found in Eq. © that 
an aperture with opening height a generates a diffraction 
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FIG. 2: (color online) a. Enhanced two-photon imaging setup 
employing different detector positions, b. By using a lens in 
the Fourier plane of the object and a two-photon absorbing 
medium our scheme is able to reconstruct the object with 
sub-Rayleigh resolution without relying on a post selection 
mechanism. 

pattern in the Fourier plane which oscillates twice as fast 
when measuring G^ 2 ^(ri,r 2 ) than when recording I(r). 
As shown above, hereby the two detectors scan the range 
in the far-field at different positions, r 2;c = —T\ m , so that 
the joint detection measurements can be performed using 
ordinary photon detectors (see Fig. [2^). However, in or- 
der to create an image of the object in the image plane of 
the lens, we could also relocate the two detectors in the 
image plane as shown in Fig. [2t>- Hereby, the condition 
T2 X = r \ x i also compatible with |r 2 J = ri x , i.e. using a 
two-photon absorbing medium, allows us to reconstruct 
the object in the image plane with sub-Rayleigh resolu- 
tion. Note that by using a two-photon absorbing medium 
the scheme does not rely on a post selection mechanism; 
in this way the scheme becomes relevant for lithographic 
applications. 

Finally, we briefly outline that our method can be ex- 
tended to different objects, e.g., to the case of an arbi- 
trary grating. Therefor, let us reconsider the expression 
of the disturbed field U (r , , ) of a single rectangular 
aperture derived in Eq. |(3} . In the case of a grating (with 
M slits of opening height a, width b and slit separation 
d), each photon may pass through either of the M slits 
before being recorded by one of the two detectors at ri 
or r 2 . For simplicity, we restrict our calculations again to 
the x-z-plane. In this case, the general expression of the 
electric field being diffracted at a grating with M slits is 
given by 

A/-1 R . M-l 

n=0 ra=Q 

1 — e jz 1 — e r 2 

= T^r- ■ ^H 10 ) 
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where, again, we made use of Fresnel and Fraunhofer 
approximations. Using this expression and by choosing 
l r 2j = r \ x ± t&S we obtain for the second order corre- 
lation function G^ (ri, r 2 , M) for a grating with an odd 
number M of slits 
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G( 2 )(r 1 ,r 2 ,M)=G^ ) (r 1 ,r 2 )- 
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cos[fc2Md^H 

' (11) 



1 - cos[fc2d^L] 



(2) 

where G (ri,r 2 ) is the second order correlation signal 
for a single slit aperture. 

The classical expression for the intensity diffraction 
pattern of a grating in case of a coherent source is well 
known 24] : it consists of a sinusoidal distribution caused 
by the grating and an envelope function Io(r) due to the 
diffraction by a single slit 



= J„(r) 
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Comparing Eq. (Ti~Tj) with (Ti"2"|) , we see that both expres- 
sions can be written as a product of an envelope function 
due to the diffraction by a single slit (G y J (i"i,r 2 ) and 
io(r), respectively) and a sinusoidally oscillating func- 
tion. However, in case of the G^ 2 ^-function the modu- 
lation frequency is twice as high as in the classical case 
corresponding again to sub-Rayleigh resolution. 

Summarizing, the proposed setup can be used to image 
a physical object with sub-Rayleigh resolution using un- 
corrected single photon sources and linear optical tools 
only. In earlier articles it was demonstrated how a 

source of single-photon emitters with a separation d < X 
can be imaged and resolved using ordinary photon de- 
tectors and joint detection techniques. In contrast, the 
scheme developed in this paper enables to resolve details 
of a distinct physical object, e.g. an rectangular aperture 
or a periodical structure, impossible to resolve by clas- 
sical far-field imaging techniques. We remark that the 
method can be implemented in various physical systems, 
e.g., with current ion trap technology [3, [13, EH- 
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